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We present a powder inelastic neutron scattering study of magnetic excitations in Ba2FeRe06, a 
member of the double perovskite family of materials which exhibit half-metallic behavior and high 
Curie temperatures. We find clear evidence of two well-defined dispersing magnetic modes in its 
low temperature ferrimagnetic state. We develop a local moment model, which incorporates the 
interaction of Fe spins with spin-orbit locked magnetic moments on Re, and show that this captures 
our experimental observations. Our study further opens up double perovskites as model systems to 
explore the interplay of strong correlations and spin-orbit coupling in 5d transition metal oxides. 
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Strong electronic correlations in the 3d and 4d tran- 
sition metal oxides (TMOs) lead to such remarkable 
phenomena as high temperature superconductivity in 
the cuprates [1 , colossal magnetoresistance in the man- 
ganites [2], and possible chiral superconductivity in the 
ruthenates [3]. In 5d-TMOs, the traditional viewpoint 
suggests that the larger spread of atomic wavefunctions 
leads to a smaller local Hubbard repulsion and a larger 
overlap between neighboring atomic orbitals, which coop- 
erate to suppress strong correlation effects. Indeed, sim- 
ple oxides like Re03 are good metals [4] . This traditional 
picture has been challenged by recent work on iridium- 
based complex oxides, which shows that the large spin- 
orbit (SO) coupling on Ir can split the X>2g crystal field 
levels, yielding a reduced bandwidth for effective j=l/2 
electrons and the re-emergence of strong correlations [5] . 
Iridates like Na2lr03, Na^raOg, Eu2lr207 and Y2lr207, 
are of great interest since they may support correlated 
SO coupled magnetism or topological phases [6HT4]. 

A distinct route to strong correlations in 5d-TMOs 
may be realized in ordered double perovskites (DPs), with 
chemical formula A2BB'06, obtained by stacking alter- 
nating ABO3 and AB'03 perovskite units. In this struc- 
ture, neighboring B(B') site ions are further apart by y/2 
(see Fig. [TJ. With 5d metal ions on the B' sites, the 
larger B'-B' distance suppresses direct 5d orbital over- 
lap, enhancing strong correlations. There is a growing 
interest in DPs such as A 2 FeRe0 6 (A=Ca,Sr,Ba) [T5ti2T] 
and Sr2CrOs06 [22j|23], with a 3d magnetic B-ion and 
a 5d B'-ion, since they exhibit high Curie temperatures. 
In addition, the half-metallic character and significant 
polarization of many DPs makes them ideal candidates 
for spintronic applications such as spin injection [24^ [25] . 
Despite this great fundamental and technological interest 




FIG. 1. Crystal structure of Ba2FeRe06- (a) Schematic of 
the crystallographic unit cell showing the relative orientation 
of Fe and Re moments, (b) Projection into ab-plane illustrat- 
ing the separation between B(B') site ions, and the exchange 
pathway J between neighboring Fe and Re sites. 

in the DP materials [26-29 , there is a significant lack of 
experimental work on their magnetic excitations. 

In this Letter, we address this issue by using inelastic 
neutron scattering (INS) on polycrystalline Ba2FeRe06 
to study the magnetic excitations in its ferrimagnetic 
state, complemented by a theoretical modelling of the 
observed spectrum. Our main results, which are sum- 
marized in Fig. [5J are as follows, (i) We provide experi- 
mental evidence of two dispersive magnetic modes in the 
magnetic excitation spectrum, showing that Fe and Re 
electrons both exhibit strong correlations and contribute 
to the magnetization dynamics, (ii) We find evidence of 
nearly gapless magnetic excitations in the inelastic spec- 
trum, indicating a weak locking of Re-moments to the lat- 
tice in the ferrimagnetic state, (iii) We discuss a minimal 
local moment model of strongly coupled spin and orbital 
degrees of freedom on Re interacting with spins on Fe, 
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which captures our experimental observations, (iv) We 
combine our results with published magnetization and X- 
ray magnetic circular dichroism (XMCD) data to obtain 
estimates of the Re and Fe moments and the effective 
Re-Fe exchange interaction. Our work further opens up 
5d-based DPs as model systems to study the interplay of 
spin orbit coupling and strong electronic correlations. 
Experiments: Total 8.6g of polycrystalline Ba2FeRe06 
sample was synthesized using the standard solid-state 
method reported previously [20j [30] . In some DPs, anti- 
site disorder (mixing of B and B' site atoms) is significant, 
and suppresses saturated magnetic moments. However, 
for Ba2FeRe06 , a large difference (~8%) in the ionic 
radii of Fe 3+ and Re 5+ seems to mitigate this problem. 
From the structural refinement of x-ray powder diffrac- 
tion data, we infer an anti-site disorder of <1%, consis- 
tent with that reported by Winkler et al. [31]. 

Neutron scattering measurements were carried out 
on the fine resolution Fermi-chopper spectrometer SE- 
QUOIA at the Spallation Neutron Source (SNS) at Oak 
Ridge National Laboratory (ORNL). Measurements were 
performed with Fermi chopper 1 rotating at a frequency 
of 300 Hz and phased for incident energies of 27 and 
120 meV. A TO chopper rotating at 180 Hz was used to 
eliminate a fast neutron background. The sample was 
sealed in an Al can and mounted on a closed cycle cryo- 
stat. Data were also collected for an empty Al sample 
can at T = 34 K, with an identical instrumental config- 
uration. The absorption corrected empty can intensities 
were subtracted from the raw data at T = 34 K to remove 
scattering from the sample environment. 

Throughout this article we use pseudo-cubic notation 
a = b = c ^ 4.01 A and index the momentum transfer Q 
in units of 1/a to aid comparison with theoretical calcula- 
tions. In our magnetic model ferrimagnetism arises from 
G-type antiferromagnetic arrangement of inequivalent Fe 
and Re moments so that magnetic Bragg peaks occur at 
the antiferromagnetic wavevector Qaf = (tt, 7r, tt) and the 
ferromagnetic wavevector Qfm = (27r, 0, 0). 

Maps of the inelastic neutron scattering intensity for 
27 meV and 120 meV incident neutron energies are shown 
in Fig. [2] (a) and (b) respectively. An inelastic feature 
emanating from Q = 1.35 A" 1 corresponding to Qaf 
is clearly resolved. The inelastic feature extends into 
two-bands of excitations with maximum intensities near 
25 meV and 39 meV. The scattering is strongest at low-Q 
and decays rapidly for increasing Q as is expected gener- 
ally from the form factor dependence for magnetic scat- 
tering. Results from our theoretical model are shown in 
Fig. [2] (c)-(d) with the best-fit parameters. 

The temperature and energy dependence of puta- 
tive magnetic scattering in Ba2FeRe06 is presented in 
Fig. [3] Bragg peaks at Q = 1.35 A" 1 corresponding 
to Qaf, and Q = 1.56 A -1 corresponding to Qfm are 
shown in Fig. [3] (a). The elastic magnetic intensity de- 
creases upon warming and the antiferromagnetic Bragg 
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FIG. 2. Comparison of spin wave model and measured pow- 
der averaged magnetic scattering intensity for T = 34 K. An 
empty Al-Can background signal has been subtracted from 
the data, (a) and (b) are neutron scattering data for incident 
energies of 27 and 120 meV respectively, (c)-(d) Powder aver- 
aged dynamic structure factor calculated from the spin wave 
model with 6j7eff T — 39 meV, and §JqrR, = 25 meV. Solid 
white lines in (a) and (b) show the dispersion relation from the 
spin wave model along (0, 0, 0)-(7r, tt, 7r)-(37r/2, 37r/2, 37r/2). 



peak vanishes above 300 K, consistent with the reported 
T c « 304 K for Ba 2 FeRe0 6 [H]. Constant momen- 
tum transfer cuts detailing the inelastic scattering emerg- 
ing from the magnetic zone center are shown in Fig. [3] 
(b). The fluctuation-dissipation theorem S(Q,E) = 
(n(E,T) + l)x"(Q,E) relates the imaginary component 
of the dynamic susceptibility x"{Q^) to the dynamic 
structure factor measured directly by neutron scatter- 
ing where n(E,T) is the Bose thermal occupation factor. 
Correcting the INS intensity by the Bose factor allows 
for comparison of the inelastic scattering across the en- 
tire 400 K temperature range on a single intensity scale. 
Two strong inelastic features are visible near 25 meV 
and 39 meV which decrease in intensity upon increasing 
temperature. The temperature, momentum, and energy 
dependence of the low-Q inelastic scattering is entirely 
consistent with expectations for scattering from powder 
averaged spin- waves. Broader examination of the data 
reveals two bands of phonon scattering which partially 
obscures the magnetic signal above 3 A -1 ; however, the 
phonon and magnetic scattering are well resolved since 
the magnetic form factor rapidly attenuates the magnetic 
intensity with increasing Q while the phonon scatter- 
ing intensity increases with Q (see Supplemental Mate- 
rial). Constant energy cuts across the low energy mag- 
netic scattering are shown in Fig. [3](c). An inelastic fea- 
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FIG. 3. Representative cuts through S(Q,E). (a) Temper- 
ature dependence of magnetic Bragg peaks at {7r,7r,7r} and 
{27T, 0, 0} determined by integrating the Ei — 120 meV data 
over the elastic line resolution [—4 < E < 4 meV]. (b) Temper- 
ature dependence of Bose factor corrected inelastic scattering 
near the magnetic zone center, (c) Constant energy cuts from 
the Ei — 21 meV data at 34 K, each energy cut was integrated 
over ±1 meV, the 7 meV cut is offset for clarity, (d) Bose fac- 
tor corrected inelastic scattering at 34 K integrated over a 
magnetic Brillouin zone representing the magnetic density of 
states, the solid line is an equivalent cut from the spin wave 
model calculation. An empty can background has been sub- 
tracted from data in (c) and (d). 



ture emerging from the antiferromagnetic zone center is 
clearly resolved within our experimental resolution down 
to at least 3 meV. The scattering intensity is strongest 
near the antiferromagnetic wavevector at Q = 1.35 A" 1 — 
where the structure factor for magnetic scattering is max- 
imized — and is small near the nuclear Bragg peak. This 
Q-dependence identifies the low energy inelastic scatter- 
ing as magnetic in origin and places an upper bound of 
3 meV for any gap in the spin wave dispersion. The 
Q-integrated inelastic intensity is peaked at the magnon 
zone boundary energy where the density of states for spin 
waves is maximized, enabling a precise determination of 
the zone-boundary energies from the powder averaged 
spectrum. The dynamic susceptibility integrated over 
the magnetic Brillouin zone is shown in Fig. [3] (d) ; scat- 
tering is strongly peaked at 25 meV and 39 meV. An 
equivalent cut from the powder averaged spin wave the- 
ory using the same parameters as in Fig. [2] (c) and (d) is 
also shown in the figure. 



Local moment model: The well-defined magnetic 
modes in Fig. [2| and the fact that the closely related 
material Ca2FeRe06 is an insulator, suggests that strong 
electronic correlations are important in Ba2FeRe06- A 
local moment model thus provides a useful vantage point 
to describe its magnetic excitations. The simplest such 
Hamiltonian is H=J Y^irr') S r -T r '—\ ^2 reRe L r -S r , with a 
nearest neighbor antiferromagnetic exchange interaction 
between the Fe spin T and the Re spin S induced by in- 
tersite tunneling. In addition, we include SOC between 
orbital (L) and spin (S) angular momentum on Re. We 
ignore SOC on Fe. This model should be broadly ap- 
plicable to many DPs with an orbitally nondegenerate 
magnetic B-site (e.g., Fe), and a magnetic B'-site with 
active t 2g orbitals (e.g., Re). 

On the Fe sites, a nominal valence assignment of Fe 3+ 
together with a strong Hund's coupling leads to a spin 
T ' = 5/2. On the Re sites, a nominal valence assignment 
of Re 5+ (5d 2 ) leads to two electrons in the t 2g orbital. 
Thus, in contrast to the iridates, not only SOC but also 
Hund's coupling (Jh) is important in determining the 
magnetic state on Re [32 j. The interaction Hamiltonian 
together with the SOC, when projected to the t 2g orbital 
[33] , takes the form 

H Re = -2 J H S 2 - J -fP - \{h • Si + t 2 • ? 2 ), (1) 

where S = si + s 2 and L = l\ + £ 2 . As seen in Fig. [^a), 
i^Re supports a 5-fold degenerate ground state over a 
wide range of J/j/A (see Supplemental Material). For 
Jh/^ ^ lj FigQb) shows that this ground state manifold 
may be viewed as made up of L = 1 and 5=1 moments 
locked into a state with total angular momentum 7? = 
L + S, with 11 = 2. 

For Ba2FeRe06 , where A ^> J , the local moment 
Hamiltonian simplifies to H e ^ = J e R ^(rr') 'R-r'Fr' , yield- 
ing an effective Heisenberg model with moments 7£, T 
on the Re and Fe sites, respectively. We find that 

J ee = J ^ + l)+5(5 ++ l)-L(L + l) . for L=s=h and 7e = 2) 

we obtain J7 e ff = J /2. We expect that the metallic na- 
ture of Ba2FeRe06 , and the concomitant carrier der- 
ealization, will lead to a smaller effective value of J 7 , 1Z 
compared to this highly localized viewpoint. 
Spin wave dispersion: The model has a ferrimag- 
netic ground state, with antiferromagnetic order of T and 
7? leading to a net ordered moment. A spin wave calcu- 
lation [34] around this ground state yields two modes, 

with energies fi±(q) = 



<S+7o — -T^Tq i <S_7o where 
S± = (J r ±1Z)/2 1 and 7 q = 2J e ^(cosq x +cosq y + cosg 2 ), 
with 70 = 6j7eff- At the ordering wavevector, Q_ has a 
gapless quadratic dispersion, while has a gap 2<S_7o. 
At T = 0, we find the dynamic structure factor for trans- 
verse spin fluctuations 



S ± (q, u) = 2tt V (G q -*S-)5(l>- n a (q) ) 



(2) 



a=± 
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FIG. 4. (a) Spectrum of H^e (in units of A) versus J#/A, with 
degeneracies indicated in brackets. For Jh = 0, the eigenstates 
at three distinct energies correspond to different ways of fill- 
ing two electrons into SO coupled single-particle states with 
angular momentum j — 3/2 and j — 1/2. When A = 0, we 
find two-particle angular momentum eigenstates 3 P, 1 D, and 
1 S in increasing order of energy. A weak SO coupling, with 
Ji//A> 1, splits the lowest 3 P manifold into 3 P 2 , 3 Pi, 3 Po- 
(b) Interaction dependence of total orbital and spin angular 
momentum, and magnitude of SO energy. For Jh > A, the 
ground state is composed of L — 1 , S — 1 moments which lock 
to yield a total angular momentum 1Z — L-\-S, with 1Z — 2. 



where G q = (70^+ — 7 q ^7£)/y <S+7q — TIZ^. As shown 

in Fig. [2|a-d) and Fig. |3jd), setting = 39 meV and 
7o7£ = 25meV in the theoretical plots leads to a broad 
agreement between the experimental data and the pow- 
der averaged theoretical result for S(q, cj), both in the 
existence and dispersion of the two magnetic modes, and 
in the (near) gaplessness of the lower energy mode. 
Spin-orbital locking on Re: For momenta with 7 q = 0, 
the spin wave dispersion yields = T ' /1Z. Since 

these momenta dominate the magnon density of states, 
we can use the ratio of the observed peak positions in 
Fig. [3] (39 meV, 25 meV), to deduce that T/Knl.G. If we 
assume that the Re moments have a pure spin origin, we 
have to set 1Z < 1. This assumption, however, yields a Fe 
spin T < 1.6, which is anomalously low — first principles 
calculations [20 , 35], a naive valence assignment of Fe 3+ , 
and the measured large saturation magnetization [30 , 
all point to a much larger Fe moment. Our observations 
thus strongly suggest that we must have 1Z > 1, indicat- 
ing a nonzero orbital contribution to the Re moment, in 
qualitative agreement with XMCD measurements. 

In order to obtain estimates of the moment sizes and 
the exchange coupling, we combine our INS results with 
previous XMCD and magnetization data. XMCD mea- 
surements indicate a significant static orbital contribu- 
tion to the magnetization on Re, with MR e b //iR e m ~ — 0.3. 
This allows us to set L~0.65', which yields S ~ 0.637^ and 
L w 0.377^. High field magnetization measurements on 
Ba2FeRe06 indicate a saturation magnetization m sat ~ 



3/i£. Together with our neutron data, this constrains 
the moment sizes to be 1Z « 1.3 and T « 2.1, yielding 
an estimated exchange coupling J e ^ ^3.1 meV. We have 
checked that including a small direct Re-Re Heisenberg 
exchange ~ 0.1^7eff slightly modifies the spin wave disper- 
sion but does not significantly affect our estimate of 7Z. 
(A large Re-Re exchange coupling leads to a dispersion 
which is not consistent with our data.) Thus, while pre- 
vious XMCD measurements on Ba2FeRe06 have shown 
that there is a static orbital contribution to the ordered 
magnetic moment on Re in the ferrimagnetic state, our 
work shows that such SO locked moments on Re also play 
a role in the low energy magnetic excitations. 
Magnetic transition temperature: We use the above 
values of the moment sizes and exchange couplings to 
estimate the magnetic T c . The nearest neighbor classi- 
cal Heisenberg model on a three-dimensional cubic lat- 
tice, with moments J 7 , 1Z on the two sublattices, has a 
mean field transition temperature IJ^TTZ. Assuming 
a quantum renormalized T c « ZJqr \J T1Z ( \J T1Z + 1), 
we estimate T c w 315K, in rough agreement with the 
measured T c expt « 304K. If one takes the limit of fully 
localized moments, setting T = 2.5 and 1Z = 2, one ob- 
tains T c w 520 K, remarkably close to that of the insu- 
lating compound Ca2FeRe06 . T c calculations retaining 
the itinerant Re electrons will be reported elsewhere. 
Structural transition and absence of spin gap: 
Ba2FeRe06 has a weak tetragonal distortion, with c/a< 
1, which onsets at the magnetic T c [19]. Since a Jahn- 
Teller distortion would lead to c/a > 1, not necessarily co- 
incident with T c , we ascribe this distortion to SOC. Going 
beyond i^ e ff, we expect a term — e X^r^r^+^r^+^r,;?) 5 
arising from the cubic anisotropy, which locks the Re mo- 
ment (and thus also the Fe spins) to the crystal axes. 
Such a magnetostructural locking term with e > ex- 
plains the observed tetragonal distortion at T c as arising 
from weak orbital order, and would lead to a spin gap 
of order e. This locking is expected to be small; on ex- 
perimental grounds since we find no clear evidence of a 
spin gap, and on theoretical grounds since it arises from 
a spin-orbit induced weak mixing of well-separated ti g 



and e g crystal field levels 



A small magnetostruc- 



tural locking term is consistent with the measured weak 
coercive field ~0.2 Tesla. 

Summary: We have used inelastic neutron scattering 
and theoretical modelling to study the magnetic exci- 
tations in Ba2FeRe06, inferring the presence of strong 
correlations and spin orbit coupled moments on Re, and 
obtaining a broad understanding of the phenomenology 
in its ferrimagnetic state. Further efforts are necessary 
to synthesize single crystals or good quality thin films 
of Ba2FeRe06 and other DPs. In future work, we will 
extend our experiments to other DP materials, and in- 
corporate the itinerant character of Re electrons in our 
theoretical modelling, both of which would lead to a bet- 
ter understanding of novel 5d-based TMOs. 
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INTERACTION EFFECTS: ATOMIC LIMIT 

For the , , d^ configuration of electrons in the 
tig orbital, we have to consider matrix elements of the 
Coulomb interaction on the same footing as the spin orbit 
coupling. The interaction Hamiltonian projected to the 
tig orbitals is given by [36 

2 



a<(3 



After some algebra, this can be reexpressed in terms of 
rotationally invariant operators as 



U-3J H 2 
-Hint = z n, 



2 ■ "tot — jj ~ tot 2 

where we assume the normal ordered form of these op- 
erators. For a d^ configuration, n tot = 2. Including 
the spin orbit coupling term leads to the effective atomic 
Hamiltonian for Re 



2JnS i 



tot 



Jh ?2 



J tot 



(4) 



#Re 



-2J H S 2 



\(£i ■ S! + £ 2 ■ s 2 ) (5) 



where L = £\ + £2 and S = si + S2 . To diagonalize this 
Hamiltonian for a d^ configuration, we write the full 
Hamiltonian in the basis \L, rri£, 5, m 8 ) corresponding to 
total orbital and total spin angular momentum. Since 
the individual orbital angular momenta £\ = £2 = 1 and 
individual spin angular momenta are s± = S2 = 1/2, we use 
a shorthand for the Clebsch-Gordan coefficients, defining 
them via 



\L,m£,S,m 8 ) \L,m £ ) <g> \S,m 8 ) 

|L,m<>= C^ a |ml,m2> 

mi ,7Tl2 

|S,m.) = £c*£|aM2> 

Sl,S2 

in terms of which the full Hamiltonian becomes 
{L' ,m^S' ,m' s \H^\L,m^S,m s ) = H L T ' ^' ° 
where 



L,m£,S,m s 



(6) 
(7) 

(8) 



(9) 



if 



L,rri£ ,S,m s 



3l,L' Ss,S' &m' g ,m £ &m' s ,m s ^L,5 
u mi,m 2 U Si,s 2 ^"*l*l^rai,ra 2 ^Si,S: 



mi,77l2,Si,S2 

yL' ,m' f 



yS'm' 



v zt/ mi + l,m 2 u si-l,S2 



v ^ u mi-l,m 2 u si+l,s 2 



cfi^L) (10) 



and 



-2J ff 5(5+l)--^L(L + l) 



(11) 



Here, we must restrict ourselves to totally antisymmet- 
ric electronic states; (L, S) = (0, 0), (1, 1), (2, 0) yield the 
allowed 15 basis states. 



When Jh = 0, we find eigenstates at three distinct 
energies corresponding to filling two electrons into single- 
particle states corresponding to a low energy j = 3/2 
manifold or a higher energy j = 1/2 doublet. On the 
other hand, when A = 0, we find i/R e has, in increasing 
order of energy, total angular momentum eigenstates 3 P, 
L L>, and L S. 

The numerically computed spectrum of H^e is shown 
in Fig. 4 of the paper. Over a wide range of Jh / A, we find 
a 5-fold degenerate ground state when spin-orbit cou- 
pling competes with Jh- For J#/A ^> 1, we can show 
that the 3 P ground states at A = split into spin-orbit 
coupled states which may be labelled by total angular 
momentum L + S = 2,1,0 in increasing order of energy, 
(corresponding to 3 P 2 3 P\ 3 Pq states with degeneracies 
5,3, 1). This shows that a local —XL-S, with L = S = 1, 
is a good description of the lowest energy manifold of 
states when Jij/A > 1. However, when Jh < A, this se- 
quence changes to 5, 1, 3 (in ascending order) suggesting 
that such a simple description fails. 



PHONON BACKGROUND 

The measured scattering intensity consists of a number 
of components including coherent nuclear and magnetic 
scattering, as well as incoherent processes. Additional 
background scattering originating from the sample en- 
vironment, and detector dark current is eliminated by 
subtracting the signal measured for an empty Al sam- 
ple can using identical instrumental configuration. The 
signal of interest is coherent scattering from magnons, 
which has a momentum dependent intensity dominated 
by the magnetic form factor. In general, the magnetic 
form- factor rapidly decays as a function of Q, thus the 
magnetic INS intensity will decrease with increasing Q. In 
contrast, both coherent scattering from phonons and in- 
coherent nuclear scattering intensities increase quadrat- 
ically with Q in a powder averaged measurement [37] . 
Any periodic modulations of the coherent phonon scat- 
tering arising from the structure factor should also in- 
crease in intensity with Q. 

A map of the inelastic neutron scattering at 34 K, for 
120 meV incident energy is shown in Fig. [5] (a). There 
are three bands of inelastic scattering, around 25, 40, 
and 70 meV which increase in intensity with increasing 
Q. We associate each of these with three phonon bands. 
The magnetic signal emerges from the antiferromagnetic 
zone center at Q = 1.35 A" 1 and extends into two bands 
with maximum intensities at 25 meV and 39 meV. 

To highlight the momentum dependence of the scatter- 
ing intensities constant energy cuts through each band 
of inelastic scattering are shown in Fig. [5] (b) - (d). 
Around 70 meV, [Fig. [5](b)], the scattering is dominated 
by phonons, here the momentum dependence of scatter- 
ing intensity is entirely described by the quadratic form 
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FIG. 5. (a) Neutron scattering intensity at 34 K for and incident energy of 120 meV. An empty Al can background has been 
subtracted from the data, (b) — (d) Constant energy cuts across bands of inelastic scattering at 70, 40 and 25 meV respectively. 
Solid black lines are a fit to I(Q) — A(E) + BQ 2 delimiting the Q-dependent contribution of phonon scattering at each energy. 



I(Q) = A(E) + £>Q 2 , where A is a constant function of 
Q parameterizing background originating from the small 
multiple scattering contributions to the inelastic scatter- 
ing. On average A is a decaying function of energy. In 
Fig. [5] (c) and (d) the overall intensity increases with in- 
creasing Q at high Q, and above Q = 3 A" 1 the scattering 
is dominated by phonons, as can be seen from the fits to 
I(Q) (solid black lines). However, below Q = 3 A -1 the 
INS intensity clearly increases above the phonon back- 
ground with decreasing Q. Furthermore in Fig. [5](d) the 
low Q scattering intensity modulation is consistent with 



the magnetic Brillouin zone. Thus, the magnetic scat- 
tering is well separated in Q from the phonon scattering, 
and the magnetic scattering is clearly identified through 
momentum, and temperature dependencies (see Fig. 3 
of the main text). We note that the two lower phonon 
modes, which are common to many perovskite materials, 
are at energies which are not far from the zone-boundary 
magnon mode energies. Further single crystal inelastic 
neutron scattering measurements are required to deter- 
mine whether this is a mere coincidence, or a result of 
magnon-phonon coupling in this material. 



